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■ Abstract 

\o : 

, Slowly rotating perfect fluid balls with regular center and asymp- 

"q | totically flat exterior are considered to second order in the rotation 

O" 1 parameter. The necessary condition for being Petrov type D is given 

for general perfect fluid matter. As a special case, fluids with a linear 
equation of state are considered. Using a power series expansion at the 
regular center, it is shown that the Petrov D condition is inconsistent 
^ . with the linear equation of state assumption. 



1 Slowly rotating fluid balls 

Rotating perfect fluid exact solutions of the Einstein equations would give 
great insights into the structure of astrophysical objects such as neutron 
stars and super-massive stars. Because of the complexity of the problem the 
usual method in trying to find such spacetimes is to make some a priory 
assumption, i. e. use an ansatz, and hope that solutions with that property 
will exist. Since the obtained solutions are very often unphysical, it would 
be valuable to have a simpler method which could decide about the physical 
suitability of the ansatz before making the effort to solve the problem in its 



1 



generality. The technique we propose in order to achieve this goal is to check 
whether the ansatz can remain valid in the slow rotation limit. 

According to the work of Hartlejl], the metric of a slowly rotating fluid 
ball can be written in the form 



ds 2 
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Here the functions A and B depend only on the radial coordinate r, determin- 
ing the spherically symmetric basis solution. It follows from the existence 
of an asymptotically flat exterior vacuum region that the first order small 
function uj depends only on the r coordinate. The asymptotic flatness also 
implies that the second order small functions h, m and k have the following 
angular dependence 

h = h + h 2 P 2 (cos$) , 

m = mo + m-iP-i (cost?) , (2) 
k = fc 2 -P2(cosi?) , 

where ho, h 2 , mo, m 2 and k 2 are functions of r and P 2 is the Legendre 
polynomial of order two. We assume that the fluid is rotating circularly and 
rigidly, and use a comoving coordinate system where the fluid velocity vector 
is proportional to d/dt. 

To quadratic order in the angular velocity, the nonvanishing components 
of a comoving tetrad can be chosen for the metric ([T|) as 

1 
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We denote the tetrad components of the Einstein tensor by G a b- The 
(0, 0) and (1, 1) components give the energy density and pressure 



V = -5 
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The only field equation is the pressure isotropy condition G u = G 22 , 

/ d 2 [A \ dB d(r 2 A) 

This equation together with an equation of state /i = and a central 
density value \i = /io determine uniquely the functions A and B describing 
some spherically symmetric perfect fluid solution [2] [3J. We note that the 
square roots are introduced only in order to combine terms together. After 
expanding the derivatives no square roots remain in ([6]) and in the rest of 
the equations. 

Calculating the field equations to first order in the small angular velocity, 
one obtains only one independent relation. The (0, 3) component of Einstein 
equation gives a second order ordinary differential equation for w, 

d ( 4 [Bduj\ „ d ( [B\ 

The only case when the rotation function u is known in terms of ele- 
mentary functions is the rotating Whittaker space-time. Among the rotat- 
ing states of the Whittaker fluid, there is the exactly known Petrov type 
D Wahlquist solution. The rotating Whittaker fluid and its matching to a 
vacuum exterior is dealt with in [1] and [5]. 

Calculating the field equations to second order in the rotational parameter 
we get differential equations for the functions h 2 , rri2, k 2 , and rriQ. The 
(1, 2) component of Einstein equations gives 

2rA-^- (h 2 + k 2 )+r (h 2 - m 2 ) — - 2A (h 2 + m 2 ) = . (8) 
The isotropic pressure condition G 22 = G33 yields, 

6A (h 2 + m 2 ) - r*B ' + 2rWA± (|) = . (9) 

After substracting 3rAB times the derivative of (jSJ) and using it again to 
eliminate the derivative of h 2 , the P 2 (cos$) part of the pressure isotropy 
condition Gu = G 22 gives 

dA dk 2 . dB dA , , 

r 2 B— - rA—h 2 - rB—m 2 + 2ABh 2 
dr dr dr dr 

+A(m 2 -4k 2 -5h 2 )-K 4 B^) =0. (10) 
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Equations (|S]), and (ITU|) together with the regularity condition at the 
center can be used to solve for the functions h 2 , m 2 and k 2 determining the 
shape of the fluid ball. The $ independent part of Gn = G22 yields the only 
field equation for ho and mo 
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Since 

dfidp = dpdfx 
d-d dr d$ dr 

is satisfied identically after applying (jSJ), (|2J) and ffTUl) . the existence of a 
barotropic equation of state // = //(p) does not bring any new condition. Of 
course, this holds only in the slow rotational limit, up to second order in 
the rotational parameter. Let us write /1 = /i + \i% and p = p + p 2 , where 
/i and p denote the value of the density and pressure in the nonrotating 
limit and \l 2 and p 2 are the second order small changes for slow rotation. 
If we assume that the equation of state remains the same for the rotating 
configuration, then /i 2 = ^£>2, and since /i 2 and p 2 is already second order 
small /i2^7 = P2^-- Calculating and writing out in detail 
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d 2 



AA(B - 1) - r 2^±^L _ Ar 2 A^B [ —V B 
dr dr \ dr 2 



(13) 



m = . 



Equations (TTTTl and (TT3|> can be used to determine the $ independent parts 
/to and mo of the second order small perturbation functions. 
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2 Petrov D conditions 



Since apart from the conformally flat interior Schwarzschild solution all static 
spherically symmetric spacetimes are Petrov type D, a natural ansatz in 
trying to find rotating exact solutions is to assume that they have some 
special Petrov type. In the rest of the paper we investigate whether this 
assumption can remain valid in the slow rotation limit. 

The electric and magnetic curvature components are defined in terms of 
the tetrad components of the Weyl tensor as follows [6J: 

Ei = C1010 , E 2 = C2020 ? -^3 — C1020 j 

n l — W010 j n 2 — <^2020 ) n 3 — <-^1020 • 

If we do not consider cosmological solutions with n = —p constant the Petrov 
type of the spacetime is D if and only if the following two conditions hold [7J 

2E\ + 5£i£ 2 + 2E 2 2 - El - 1H\ - bH x H 2 - 2E\ + E\ = , (14) 
4£?ifTi + 5£i# 2 + 5E 2j ffi + 4E 2 H 2 - 2E 3 H 3 = 0. (15) 

Equation (|T5l) is always identically satisfied to second order, while ( JT4")) gives 
the simple condition 

{hl - m2) = rV(±±) 2 . (16) 

This condition is derived without assuming any specific equation of state for 
the matter. In order the spacetime to be Petrov D, equation ( |T6i) must be 
consistent with (|6l), 0, (j8]h (|9|) and (fTU|). 
If 

dB , „. , _. 

r— + 2(1-5) =0 (17) 

then B = 1 — cir 2 , where c\ is some constant. Substituting into ([6]), one gets 
a second order differential equation for A, the general solution of which is 

A = c 3 (Vl-cir 2 + c 2 ) 2 . (18) 

This A and B describes exactly the interior Schwarzschild solution, for which 
it has been shown that cannot remain Petrov type D when set into slow 
rotation [7]. This can also be seen easily in the present formalism in the 
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following way. From (1161) it follows that uo = c^A/B, where C4 is some 
constant. However, this uo is not a solution of the first order field equation 

It follows that for all rotating Petrov D spacetimes + 2(1 — B) is 
nonzero. In this case equations ffT6l) and (Q can be solved algebraically for 
h 2 and m 2 , and then a linear combination of (jSj) and (fTUj) can be solved 
algebraically for fc 2 . Substituting back to (jSj) again, one obtains a lengthy 
differential equation containing only A, B, u and their derivatives. Since this 
equation is 293 terms long, we will refer to it as 6293 from now on. 

It is possible to eliminate u from 6293 and its r coordinate derivatives. One 
first obtains a lengthy equation linear in ^ by taking a linear combination of 
e 2 93 and its first derivative. It is also possible to calculate a similarly linear 
equation by combining 6293 and its second derivative. A linear combination 
of these two equations gives an enormous differential equation which does 
not contain to at all. Collecting all terms to one side, this equation turns 
out to be the product of a 65 terms long and a 2446 terms long polynomial, 
which we call e§§ and e2446 respectively. These equations are polynomial in 
A, B and the derivatives of B. The shorter equation contains up to third 
derivatives of B, while the longer contains even fifth derivatives. It seems no 
other way to avoid these higher derivatives than to use some specific equation 
of state for the fluid. 

It is also possible to completely eliminate u by taking combinations of 
6293 and only its first derivative. The resulting equation, however, turns out 
to be even longer, the product of and a 6774 terms long polynomial, e6774- 
There does not seem to be any easy way to combine e6774 with e2446 in order 
to get a shorter equation. 

The field equation (jBj) and e 65 or e 2 446 could be, in principle, solved for 
the functions A and B. This indicates that or e2446 essentially plays the 
role of an equation of state for the fluid. Because of the complicated form of 
these polynomials its seems hopeless to use (jlj) and (J5j) to get an equation of 
state in the form \x = A possible way to proceed is to check wether 
or e 2 446 may be consistent with some simpler but physical equation of states. 

3 Linear equation of state 

In the rest of the paper we investigate slowly rotating perfect fluid balls with 
the equation of state \i — ap + {3, where a and (3 are some constants. Using 
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(jlj) and (j5J), we can write the equation of state as 



1 d(rB) BdirA) „ 2 



We will use this equation to eliminate the first derivative of the function B. 

dA 
dr 



After separating two factors, /i = 4n and 



/ 2 =( a + l)£^ + A(/?r 2 -a-l) , (20) 

6293 becomes 364 terms long. 

If /i = then A is constant, and the general solution of (JSl) is £> = 
1 — cr 2 , where c is some constant. This is the Einstein static universe with 
\i = —3p = 3c. Since in this paper we are only interested in spacetimes 
describing isolated rotating bodies with regular center and a p = surface, 
we do not study the rotating version of this solution. 

The /2 = case also contains only cosmological solutions. If a+1 nonzero 
and rA is not constant then we can use ( 1201) to express B. Substituting into 
the field equation (jHJ) we get 

(a + 1 - /3r 3 ) ^ + 2(a + l)f + 2/3rA = . (21) 
The general solution of this equation is 



A= C -± + c 2 

r 



1 l 2 

1 r 

3(a + l) 



(22) 



where C\ and C2 are constants. The energy density and pressure, using 
and ([5]) , turns out to be constant 



VL = -p = S— . (23) 
a + 1 

However, this spacetime is regular at the center r = only if c\ — 0, in 
which case it is just the de Sitter solution. The rA constant subcase is 
always irregular at the center. If we look at the a = — 1 subcase, then from 
f'2 = it follows that (3 = 0, and the equation of state fTl9|) shows that B/A 
must be constant. Substituting into the field equation, it turns out that the 
only solution with a regular center is the de Sitter solution. Since we are not 
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investigating cosmological solutions in the present paper, in the following we 
assume that f% and / 2 are nonzero. 

Substituting for the derivative of B from (fT9l) and dropping an /i/ 2 fac- 
tor, becomes 80 terms long, but after dividing by (a + 3) 2 f 1 3 f 2 3 , the 
polynomial e2446 reduces to 574 terms. The a = — 3 case is the Whittaker 
solution. One of its rotating generalizations is the Petrov D Wahlquist so- 
lution. All other rotating versions of it belong to the algebraically general 
class. Unfortunately, the 574 term equation is still too long to prevent us 
from eliminating A or B by taking derivatives and linear combinations of 
them. 

A particularly simple case is the equation of state /i = p with a = 1 and 
(3 — 0. Then e 293 is only 36 terms long, e 65 becomes 8 terms, and e 2 446 is 
29 terms long. Taking repeated linear combinations of the equation e 2 446 (or 
e^) and its first derivative in order to eliminate the highest power of 44, in 
the end one gets a product of long polynomials containing only B. Since the 
coordinate r does not appear, this means that B must be constant. Solving 
the equation of state \i = p for 4^ and substituting into the field equation §6§ 
it turns out that B can only be 1 or |. The B = 1 case is just the Minkowski 
spacetime. For B = \ one obtains A = cir 2 , where c\ is some constant. 
However, since /x = p = ^2, this solution is singular at the center. This 
example shows how important it is to consider the regularity conditions at 
the center of the fluid ball. 

4 Expansion around a regular center 

A spherically symmetric spacetime described by the functions A and B is 
regular at the center of the fluid ball if the limit of B is 1 and the limit of A 
is a positive constant at r = 0. The B — > 1 condition ensures that the surface 
per radius ratio of small spheres is the proper flat space value. 

In order to know if there are Petrov type D rotating perfect fluids with a 
linear equation of state one has to decide whether there are solutions A and 
B of the field equation ([6]) and the linear equation of state condition ( |T9l 
which also satisfies the Petrov D consistency condition or e 2 446. Even 
though e 2 446 is only 574 terms long for the linear equation of state case, it 
appears to be an extremely difficult task to build up some shorter solvable 
combination of e 2 446 and its derivatives. A computationally feasible way to 
get around this problem is to calculate a power series expansion solution of 
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flSD and (TTU|> and check if it does satisfy the consistency condition e 65 or e 2 u6- 
Since we are interested in configurations which are regular at the center of 
the fluid ball, it is practical to do the expansion around the center r = 0. We 
take the expansion of the functions A and B in the form 

oo oo 

A = J2Ar l , B = J2B ir \ (24) 

i=0 i=0 

where Ai and Bi are constants. Regularity at the center requires that B = 1 
and also implies that Ai and Bi are zero for odd i. Then two freely specifiable 
parameter remains in the solution of the equations (jSD and (I19p . which we 
choose to be Aq and Bi- For convenience we introduce the notations Aq = a 
and B 2 = b for them. Substituting the expansion fl24l) into (jSJ) and f|T9|) . 
the rest of the expansion coefficients can be calculated recursively. To give a 
picture about their structure we list the first few: 

A 2 = ~[b(a + 3)+0\ (25) 

A * = -^[% + 3)+/3][3&(a + l)+/3](a-5) (26) 
80a z 

Ba = -^-[b(a + 3)+/3][3b(a + l)+(3} (27) 
20a 

A 6 = --^—[4b{a-l)(a-9) + P{3a-23))B 4 (28) 
168a 

B 6 = -U&(4a 2 + 9a + 15)+ / 3(3a + 5)]5 4 (29) 
28a 

As = -^Jf|^[ &2 («- 13 )(«- 3 )( 61 « 2 + 6 « + 45 ) 

+6ty3(15a 3 - 211a 2 + 321a + 195) (30) 
+p 2 (3la 2 - 354a + 195)]5 4 

B 8 = 1 2 [6 2 (61a 4 + 200a 3 + 618a 2 + 864a + 945) 

+66/?(15a 3 + 44a 2 + 96a + 105) (31) 
+/3 2 (31a 2 + 96a + 105)]5 4 

A " = -I33» |6,,3(a - 5)(a - 17)(629a4 - 350a3 
+768a 2 - 378a + 675) + 6/? 2 (8275a 5 - 173868a 4 

+618128a 3 - 44841a 2 + 111645a + 344250) (32) 

+66/3 2 (962a 4 - 18225a 3 + 46297a 2 + 28065a 
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+19125) + /3 3 (1261a 3 - 20672a 2 + 33285a 
+ 12750)] J B 4 . 

It turns out that in order to be able to decide about the Petrov D condition 
one has to calculate this expansion at least up to % = 16. The calculations 
were performed using the algebraic computer language Reduce [8]. 

Substituting the power series expansion forms of A and B into the Petrov 
D consistency condition e 65 or e 2 446, they have the form 

oo oo 

e 65 = C i T% ' e 2446 = DiT% ' ( 33 ) 
i=0 i=0 

where the constant coefficients Ci and Di can be calculated in terms of Ai 
and Bi. 

The first nonvanishing coefficients in are 

C 6 = — [b(a + 3) + (3f [36(a + 1) + (3} (34) 
a 



c 

" 6(22a 2 - 37a - 189) + 7(3(2a - 9) (35) 



°io = ^[/3 2 (144a 2 - 1075a + 1635) 



50a 
2800a^ 

+b(3(A67a 3 - 2624a 2 - 1509a - 9810) (36) 
+6 2 (355a 4 - 1348a 3 - 5430a 2 + 5148a + 14715)] . 



In order to be Petrov D in the = subcase, all Ci must vanish. Since all 
coefficients contain the factor Cq, this is equivalent to the C& = condition. 
The constant a cannot be zero, since A would also be zero then. 
If b(a + 3) + (3 = then A = a constant, and 

B = l- J?—r 2 . (37) 
a + 3 

This is the Einstein Universe with 

A* = -3p=-^r. (38) 
a + 3 



If 36(a + 1) + j3 = 0, then 



B = - = l -^zr/ ' ( 39 ) 

a 3(a + l) 
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which is the de Sitter solution with 



a + 1 

Hence we have proven that the condition is inconsistent with ([6]) and 
( TT9j) if one considers fluid balls with a regular center and asymptotically flat 
exterior. 

The first nonvanishing coefficients in e 2 446 turn out to be 
72a 8 

D 12 = ^—Jb(a + 3) + P} 2 [3b(a + l)+(3] 2 

[b(lla + 21) + 7(3} [b{5a 2 + 54a + 9) + (3(5a + 3)] (41) 
D u = -^- 5 [b(a + 3) + (3} 2 [3b(a + 1) + (3} 2 

[6 3 (5083a 5 + 59975a 4 + 76170a 3 + 181170a 2 
+443691a + 16119) + 6 2 /3(11954a 4 + 102947a 3 (42) 
+204195a 2 + 378369a + 16119) 
+6/5 2 (9071a 3 + 54535a 2 + 104349a + 5373) 
+/3 3 (2200a 2 + 9175a + 597)] . 

All coefficients Di contain the same two factors as Cj, however, as we have 
seen earlier, these factors are only vanishing for cosmological spacetimes. 
Calculating also D 16 , the non-common factors in these coefficients can be 
shown to be mutually contradictory. Hence we have proven 

Theorem: A slowly rotating perfect fluid configuration with linear equa- 
tion of state, regular center and asymptotically flat exterior cannot be Petrov 
type D. 

Our method of deciding about the Petrov type can be relatively easily 
generalized to perfect fluids with non-linear equations of states as well. The 
algebraic calculations become more complicated, but the essential point is 
that there is no need to solve differential equations in the process. Another 
important future work would be to check if, apart from the Petrov types, 
various other physical properties or generalized symmetries may or may not 
remain valid in the slow rotation limit. That research could yield the first 
step towards finding a physically realistic rotating perfect fluid exact solution. 
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